Numerical study of resistivity of model disordered three-dimensional metals 
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We calculate the zero-temperature resistivity of model 3-dimensional disordered metals described 
by tight-binding Hamiltonians. Two different mechanisms of disorder are considered: diagonal 
disorder (random on-site potentials) and off-diagonal disorder (random hopping integrals). The 
non-equilibrium Green function formalism provides a Landauer-type formula for the conductance of 
arbitrary mesoscopic systems. We use this formula to calculate the resistance of finite-size disordered 
samples of different lengths. The resistance averaged over disorder configurations is linear in sample 
length and resistivity is found from the coefficient of proportionality. Two structures are considered: 
(1) a simple cubic lattice with one s-orbital per site, (2) a simple cubic lattice with two d-orbitals. 
For small values of the disorder strength, our results agree with those obtained from the Boltzmann 
equation. Large off-diagonal disorder causes the resistivity to saturate, whereas increasing diagonal 
disorder causes the resistivity to increase faster than the Boltzmann result. The crossover toward 
localization starts when the Boltzmann mean free path £ relative to the lattice constant a has a 
value between 0.5 and 2.0 and is strongly model dependent. 



Saturation of resistivity in metallic compounds^ as well 
as its absence^ is an interesting phenomenon which is 
far from fully understood. Some compounds saturate at 
the levels predicted by the lofFe-Regel condition {£ = a, 
where a is lattice constant and £ is the mean free path), 
others saturate at much larger levels (i.e. higher resis- 
tivity), and there are some that do not saturate at all. 
This diversity in behavior received substantial attention 
from theorists. For example Millis et alw- applied dynam- 
ical mean-field theory to calculate resistivity of electrons 
coupled to phonons and static disorder. Gunnarsson et 
dlA studied several rather realistic models with different 
forms of electron-phonon coupling using quantum Monte- 
Carlo method. They observed saturation of resistivity 
in the case of phonons coupled to hopping matrix ele- 
ments. In an attempt to understand the mechanism of 
saturation we chose to study transport properties of sim- 
ple models of metals with static disorder (as opposed to 
models closely reproducing reality). 

The resistivity was calculated using the Landauer-type 
formula for the zero - temperature linear response of 
mesoscopic systems, which can be derived in the frame- 
work of non-equilibrium Green function formalism^: 

G = '^Tr[TLG'°'TRG^''] 
h 

where G"'''* and G^'^^ are retarded and advanced Green's 
functions of the system of interest, and T ^ are matri- 
ces describing the effect of contacts on the system. 

The formula above is suitable only for finite samples, 
whereas our aim is to calculate the resistivity of bulk 
disordered material as a function of strength of disorder. 
The solution is to calculate resistance of several samples 
of different lengths and then extract bulk resistivity p 
from the data using the formula R = Lp/A {L is length. 



A is cross-section). 

Let us describe the set-up of calculation for an in- 
dividual sample. We consider a sample consisting of 
Nx X Ny X Nz unit cells (for our simple cubic exam- 
ples, each unit cell contains one atom.) The sample is 
placed between two semi-infinite contacts of the same 
cross-section Ny x Nz- Both contacts and sample have 
the same crystal structure and are described by a tight- 
binding Hamiltonian with the same parameters. Then 
tight-binding parameters for atoms inside the sample are 
randomly changed from their initial values according to 
rules given later and the resistance of the disordered sam- 
ple is calculated. 

Our sample is not periodic in the direction of current 
flow X, but, in order to decrease the effect of boundaries 
on the results, periodic boundary conditions (with pe- 
riod Ny and N^) are used in the perpendicular directions 
Y and Z . Then standard fc-vector formalism applies in 
these two directions. For a given fc-vector, Hamiltoni- 
ans of the sample and contacts are constructed and the 
conductance is calculated. Then the conductance is aver- 
aged with equal weights over fc-points on a uniform grid 
(a 6 X 6 grid was used, with a new random Hamiltonian 
at each fc-vector.) 

In the limit of small disorder strength, three- 
dimensional transport should be accurately described 
by the linearized Boltzmann equation. Therefore, it is 
useful to compare our numerical results with resistiv- 
ity obtained from the Boltzmann equation. We did not 
solve this equation exactly, but used instead the stan- 
dard procedure of a displaced Fermi-Dirac distribution^ 
F{k) ~ fFDik + eEr/h), where the displacement eEr/h 
to variational accuracy is given by: 

h _ 2^ ^kk' l^fcfc'|^(«fc -Vk')'^S{ek - eF)S{ek' - £f) 
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where hvk = ^k^k, Vkk' is the matrix element of scatter- 
ing potential calculated in Born approximation and the 
bar indicates the ensemble average. 

The first case is a simple cubic crystal, with one atom 
per unit cell and one s-orbital per atom. There are two 
tight-binding parameters — the energy level eq (diagonal 
element of Hamiltonian) which is taken to be zero and 
the hopping integral t — 1 between first nearest neighbors 
(all other hopping integrals are neglected). We consider 
a half-filled band. Two types of disorder are possible in 
this model — diagonal and off-diagonal. 

In the case of the diagonal disorder, the hopping pa- 
rameter t is kept constant throughout the sample but the 
energy level e is changed randomly according to the for- 
mula: e = + ^, where ^ is a random variable distributed 
uniformly in [—W/2,W/2]. W serves as the measure of 
disorder strength. The scattering potential \ Vkk' |^ is then 
W^/12. 

In this calculation as well as in all others, the cross- 
section of samples is 9 x 9, and the lengths used are 5,6,7 
and 8. For each length, 36 configurations of disorder are 
created, the resistance of the sample for each configura- 
tion is calculated and then averaged over configurations. 
This procedure is repeated for different L and W. In 
order to find the bulk resistivity, resistance vs. L is plot- 
ted for each W and resistivity is found from the slope of 
resulting line. 

Final results arc shown on Fig 1. The resistivity pB 
calculated using linearized Boltzmann equation is plotted 
on the same graph for comparison. The resistivity de- 
pends linearly on in a quite large range of W (up to 
W = 8). At larger W it deviates upwards from the Boltz- 
mann resistivity pb- The dotted curve shows the mean 
free path £ —< >i/2 j- calculated from Boltzmann 
resistivity vs . It can be seen that deviation of resis- 
tivity from linear starts approximately when (./a 0.5. 
This result reproduces the earlier work done by Nicolic 
and AUen^ (with a different computer code.) 

In the case of off-diagonal disorder, the energy level e is 
kept constant whereas hopping elements t inside the sam- 
ple are changed randomly according to: i = 1 + ^, where 
^ is again a random variable distributed uniformly in 
\—W/2,W/2\. The calculational procedure is otherwise 
the same. Results for off-diagonal disorder are shown on 
Fig 2. It can be seen that resistivity depends linearly on 
W"^ up to = 2. At larger W, when Ija^ 1, the resis- 
tivity deviates from linear, and at = 4 it starts to sat- 
urate at some level. Therefore in the case of off-diagonal 
disorder there is no metal-insulator transition (in agree- 
ment with the statement of Antoniou and Economou'^). 
It should also be noted that our result is similar to the re- 
sults of Calandra and GunnarssonS, obtained for a more 
realistic model of disorder induced by electron-phonon 
interaction. 

We also studied the combination of diagonal and off- 
diagonal disorders, in order to see how interplay be- 
tween them would affect resistivity. Diagonal disorder 
was kept constant and off-diagonal disorder was gradu- 
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FIG. 1: Resistivity for a single s-band model with diagonal 
disorder. Solid line is the Boltzmann resistivity for the same 
model. Dashed line shows the Boltzmann mean free path vs 



CD 
CM 



Q. 

>^ 
■> 

'05 
CD 
CC 




Resistivity 
Boltzmann resistivity 
Boltzmann mean-free path 



cd 

Q. 

CD 
CD 



c 
CO 



5 10 15 

Disorder strength (W/t) 



20 25 

2 



30 



FIG. 2: Resistivity for a single s-band model with off- 
diagonal disorder. Solid line is the Boltzmann resistivity for 
the same model. Dashed line shows the mean free path vs 



ally increased. Some results are shown in Fig 3. It is clear 
that the presence of diagonal disorder can change the re- 
sistivity vs. off-diagonal disorder dependence. At the 
diagonal disorder Wdiag = If resistivity decreases with 
increasing off-diagonal disorder. At Wdiag — 6 resistivity 
displays nonmonotonic behavior. In all cases resistivity 
saturates eventually at approximately the same level as 
in the absence of static diagonal disorder (see Fig 2). 

Next we consider a simple cubic structure with one 
atom per unit cell, but two d-orbitals per atom, namely 
— y"^ > and |3z^ — > {Eg orbitals). As in the 
previous model, only first nearest neighbor hopping is 
considered. In this case, all hopping elements can be ex- 
pressed in terms of two coupling parameters (dda) and 
(ddS) which are taken to be 0.051 and 0.003eT^ respec- 
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FIG. 3: Resistivity for a single s-band model with combi- 
nation of disorders. Diagonal disorder is kept constant, off- 
diagonal disorder varies. 
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FIG. 4: Resistivity for tlie Eg model with diagonal disorder. 
Solid line is the Boltzmann resistivity for the same model. 
Dashed line shows the mean free path vs W'^. It can be seen 
that resistivity deviates from Boltzmann expression when 1 
becomes close to lattice constant a. 



tively. Both orbitals have energy levels set to zero. 

First, diagonal disorder with a uniform distribution 
was considered. The two Eg orbitals are given inde- 
pendent random diagonal energies. Disorder strength W 
is measured in terms of the full bandwidth, which for 
present choice of parameters equals 0.3eV. The resulting 
graph of resistivity vs is shown in Fig 4 along with 
the Boltzmann resistivity. As in the case of a single or- 
bital model, the resistivity is initially linear in and 



then begins to deviate upwards when ija ^ 1, indicating 
the presence of a metal-insulator transition at larger W . 

Off-diagonal disorder in the case of two d-orbital model 
can be created in numerous ways since there are several 
different hopping elements in the Hamiltonian. We chose 
disorder strength to be proportional to the hopping cle- 
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FIG. 5: Resistivity for the Eg model with off-diagonal dis- 
order. Solid line is the Boltzmann resistivity for the same 
model. Dashed line shows the mean free path vs VK^. It can 
be seen that resistivity deviates from Boltzmann expression 
when 1 becomes close to lattice constant a. 



ment itself, or more specifically: t — io = *Ci where ^ is 
defined above. Results arc shown in Fig 5 and are simi- 
lar to those of the single s-orbital model. The resistivity 
starts to deviate from Boltzmann- like when tj a 2, and 
to saturate when £ conies close to the lattice constant, in 
agreement with loffe-Regel condition. 

The main conclusion to be drawn is that the key 
factor to resistivity saturation in metals is strong off- 
diagonal disorder. The combination of diagonal and off- 
diagonal disorder can produce various types of the resis- 
tivity behaviour — decreasing, increasing or nonmono- 
tonic. Saturation is found even in the case of a sin- 
gle band, which is in disagreement with the statement 
of Allen and Chakraborty^° that multi-band structure 
is essential for saturation. The value of £/a at which 
Boltzmann theory starts to break down varies surpris- 
ingly strongly, from 0.5 to 2.0, for the four models con- 
sidered here. 
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